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Abstract: This paper considers several distinct mathematical and computational tools, 
namely complexity, dimensionality-reduction, clustering, and visualization techniques, 
for characterizing music. Digital representations of musical works of four artists are analyzed by 
means of distinct indices and visualized using the multidimensional scaling technique. The results 
are then correlated with the artists’ musical production. The patterns found in the data demonstrate 
the effectiveness of the approach for assessing the complexity of musical information. 
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1. Introduction 


The relationships between music and mathematics have been studied for long [1,2]. However, 
it seems difficult to find a single model for describing a musical work, in spite of it being recognized 
that we have a glimpse of mathematical structures underneath all types of music [3,4]. A musical work 
can be represented as a set of one or more time-sequenced digital data streams, reflecting a given time 
sampling of the original musical source. If a single (‘mono’) digital data stream is adopted, then a 
musical work is represented by a time series (TS), where each sample is a signed floating-point value. 

Complexity is one important characteristic of a TS and embeds a description of properties, such as 
chaoticity, fractality, regularity, and memory [5,6]. In other words, while various properties can 
describe specific aspects of the TS, the complexity constitutes a general quantitative estimation of 
their characteristics [6]. Therefore, complexity has become an increasingly prevalent estimator in 
analyzing TS produced by complex systems, such as in economics [7], finance [8], geo [9], life [10], 
and social [11] sciences, with the objective of finding the fundamental principles that govern the 
systems’ behavior [12]. There are no definite guiding rules to the interpretation of the complexity 
measurements. In general, low complexity indicates that the observed system is more likely to follow 
some kind of deterministic process that can be finely captured. On the other hand, high complexity 
represents some data dynamics that are more unpredictable and difficult to understand [6]. 

A variety of complexity indices has been adopted for tackling art, namely entropy [13,14], 
Kolmogorov complexity [15,16], fractal dimension [17,18], and others [19,20]. Despite that some of 
these tools may be correlated, they capture different aspects of the system and, therefore, complement 
each other [21]. Specifically for the case of music, we can mention the work of Simonton [22] who 
studied 15,618 themes of classical music and established a connection between melodic complexity 
and popularity. Eerola and North [23] analyzed the melodic complexity of Beatles’ songs and 
observed an increasing trend over time. Additionally, they noted some kind of correlation between 
complexity and the songs’ popularity. Herrera and Streich [24] explored the relationship between the 


Entropy 2019, 21, 669; doi:10.3390/e21070669 www.mdpi.com/journal/entropy 


Entropy 2019, 21, 669 2 of 18 


detrended fluctuation analysis exponent and several manually annotated semantic labels of musics. 
These researchers concluded that there was a link between the exponent and the idea of ‘danceability’. 
Li and Sleep [25] investigated the topic of melody classification using a similarity metric based on 
the Kolmogorov complexity. Ribeiro et al. [26] analyzed 10,000 songs using the complexity-entropy 
causality plane. The results indicated that this representation could not only discriminate, but also 
allow a quantitative comparison between songs. Several other works address the use of distinct 
complexity indices and the adoption of alternative approaches [27-30]. 

In this paper we take advantage of the synergies between several tools, namely complexity, 
dimensionality-reduction, clustering, and visualization techniques, for studying the musical repertories 
of three singers and one band of different nationalities and ways of being: Frank Sinatra, Rolling Stones, 
Johnny Hallyday, and Julio Iglesias. In a first phase, the original musical sources are converted into 
‘mono’ digital format and processed using eight distinct complexity indices. Moreover, the results are 
correlated with the periods of the artists’ careers. In a second phase, a multidimensional scaling (MDS) 
algorithm is adopted for visualizing complexity. The MDS processes the dissimilarity information 
calculated with the arc cosine and Canberra distances between the complexity indices, and the loci 
generated are interpreted under the light of the emerging patterns. 

Considering these ideas, Section 2 introduces the mathematical background. Section 3 analyses 
the musical portfolio of four well-known musicians in the perspective of the eight complexity indices. 
In addition, the MDS-generated maps are interpreted having in mind the evolution of the artists’ 
careers. Finally, Section 4 presents the conclusions. 


2. Mathematical Background 


2.1. Entropy 


The information theory [31,32] has been successfully adopted in the study of complex systems [21,33]. 
Let us consider a discrete random variable X with sample space {xj,...,Xj,...,Xm} and 
probability distribution P(X). The Shannon entropy, H, of X is given by: 


- EPG i) log P (x;). (1) 


The Jensen-Shannon divergence (JSD) measures the dissimilarity between two probability 
distributions P(X) and P(Y) and is defined as [34]: 


1[™ M M 
JSD[P(X) || PO] =5 Y P(xi) log P(xi) + )) P(yi) log P(yi) =a zi)log P(zi), (2) 
i=1 i=1 


where X and Y are random variables with sample spaces {x1,...,X;,-..,xm}and {y1,...,Yi,---,ym}, 
respectively, and Z = 5(X +). 


2.2. Permutation Entropy 


Different entropy formulations and entropy-based indices have been proposed for data 
characterization [35-39]. The permutation entropy (PE) was originally proposed to assess the 
complexity of TS [40]. Let us consider a TS consisting of a series of real-valued samples {xn : 
n = 1,...,N}. We define the parameters d,t € N, where the embedding dimension, d > 2, 
and the embedding delay, tT > 1, represent the length of the TS partitioning sequences and the 
separation time between their elements, respectively. Let us denote by ¥ = {II1,...,I1q} the set of 
all possible permutations of the ordinals {1,...,d}, and by [Z] the Iverson bracket [41], such that: 
r= 1, if Zistrue 

0, ifZis false © 
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The procedure for calculating PE can be outlined as follows: 
1. Foreachn =1,...,K, with K = N — (d—1)t, 


1.1 compose the sequence {Xn,%n47,--+,X sae ts 
Xn+t +++ Xn4(d-1)r]. 
ae d ’ 
1.3 sort the array by increasing order of the elements in the first row; 
1.4 denote by 71; the sequence of numbers in the second row of the sorted array; 


1.2. construct the 2 x d dimensional array 


2. Compute the probability distribution sD: where W is a random variable with sample space 
{wy ,...,Wj,...,Wa} and P(w;) ee ee |, fort = 1pss.pdl 
= 
3. Calculate PE as ’ 


The permutation entropy PE lies in the interval 0 < PE < 1. The minimum value PE = 0 
indicates that the TS is regular, or predictable, while the maximum value PE = 1 corresponds to a 
random TS. The embedding dimension must be chosen such that N > d! in order to obtain reliable 
values of PE. For practical purposes, the values d € {3,...,7} and tT = 1 are recommended [40]. 


2.3. Statistical Complexity 
Another complexity index is the statistical complexity, C, given by [42,43]: 


C =<. JSD [P(W) || P(U)] - PE, (4) 


where U is a random variable with sample space {t1,...,Uj,...,Ug)}, probability distribution P(U), 
and P(u;) = 7,80 that: 


a 
= 


k = max {JSD [P(W) || P(U)]} = log (d! + 1) + log d! — 2 log (2d!) (5) 
is a normalization constant. 

The statistical complexity, C, depends on a probability distribution associated with the system, 
P(W), and on the uniform distribution, P(U). Therefore, for a given PE, there exists a range of possible 
values of C. Indeed, the index C provides additional information not captured by the index PE, since it 
quantifies the existence of correlational structures in the data [42,44]. 


2.4. Kolmogorov Complexity 


The Kolmogorov complexity, K(V), of an object ¥ provides a measure of information that is 
independent of any probabilistic assumptions about the data sequences in . The measure K(“) 
is defined as the size of the shortest program that, given an empty object at its input, computes V 
in a universal computer and then stops [45,46]. The exact value of K(#’) is not computable [45,46]. 
Therefore, approximation schemes are used to obtain its upper bounds, such as the Lempel-Ziv [47], 
linguistic [48], and compression-based [49] methods. 

Lossless compression algorithms approximate K(4’) from the size of the compressed object, 
K(&) & size[®(¥)], where ®(-) denotes the compression algorithm [46]. However, for obtaining a 
good approximation, the compressor has to be ‘normal’, meaning that, given ¥ and the concatenation 
of X with itself, VV, the compressor must generate compressed objects such that size[®(V)] ~ 
size[P(X X)| [46]. Moreover, for obtaining a complexity index that is independent of size[’] we adopt 
the complexity ratio, CR, given by: 
size[®(1V)] 


as size|¥] 


(6) 
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2.5. Multidimensional Scaling 


Clustering and visualizing data with a large number of attributes is overly important in 
science [50-54]. The MDS is a computational technique for dimensionality-reduction, clustering, 
and visualization of multidimensional data [33,55-57]. Given a set of objects x;,i = 1,...,L,ina 
r-dimensional space, and a measure of dissimilarity between the pair i and j, 6;;, the procedure starts 
by calculating an L x L symmetric matrix, A = [6;;], of object-to-object dissimilarities. The matrix A is 
the input to the MDS computational algorithm. In fact, MDS represents objects by means of points 
located in a q-dimensional space (q < r) at distances 6;;. To accomplish this, the MDS iterates multiple 
configurations and calculates the matrix of distances © = [6;;] that minimizes a fitness function. 
A widely used fitness function is the raw stress: 


n= [y= fp], (7) 
where f(-) is a linear or non-linear transformation. 

The MDS interpretation is based on the patterns of points emerging in the MDS locus. Two similar 
(dissimilar) objects are shown as two points that are close to (far from) each other. Therefore, we can 
translate, rotate, and magnify the locus to have a good visualization, because the object-to-object 
distances remain identical. The MDS axes have neither units, nor special physical meaning. 

The MDS quality can be quantified by means of the Shepard and stress plots. The Shepard 
diagram compares 6;; and 6;;, for a particular value of q. A narrow scattering of the points represents a 
good fit between 6;; and 4;;. The stress diagram represents the locus of R versus q. Usually, we adopt 
q = 2 or gq = 3, because such values allow a direct visualization and establish a compromise between 
achieving low values of R and gq. 


2.6. Musical Sounds 


In the context of this study, a musical work is a TS, ¥ = {x, : i = 1,...,N}, representing 
the arithmetic average of two data streams that result from sampling the original musical source at 
F, = 44.1 kHz. 

Using the discrete Fourier transform we can express the TS in the frequency domain, resulting in: 


V=H{teik=laaNP= FIA}, (8) 


N 2a 
VY = y ane tN D1) (9) 


n=1 


where j = /—1 and F{-} is the Fourier operator. Usually, we represent only the first half of the 
spectrum versus frequency, f, or angular frequency, w = 27f, by considering k = 1,...,[] and 
fakg/(7l. 

The musical sounds have a strong variability, making difficult their quantitative characterization 
through a single index. Therefore, often several distinct indices are used in the time and frequency 
domains to capture the rich information embedded in the signal. Figure 1 illustrates the musical 
work ‘LA is my lady’ by Frank Sinatra using its TS and amplitude spectrum representations, ¥ and 
|Y|, respectively. 
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Figure 1. Representation of the the musical work ‘LA is my lady’ by Frank Sinatra in the time and 
frequency domains: (a) Time series (TS), 1; (b) amplitude spectrum, |)|. 


A variety of features have been proposed for characterizing musical sounds in terms of 
their dynamics, rhythm, timbre, pitch, and tonality. Music feature extraction involves many 
signal-processing techniques and forms the basis for many automatic classification algorithms. Several 
toolboxes are currently available for music and sound feature extraction, such as the MIRtoolbox [58], 
pyAudioAnalysis [59], and Librosa [60]. The toolboxes often provide not only a set of base features 
that capture various temporal, spectral, and spectrotemporal properties of the musical signal, but also 
a considerable number of descriptors derived from the base features by means of descriptive statistics. 
Typically, all toolboxes provide onset detection, pitch tracking, mel frequency cepstral coefficients 
(MFCC), chroma, and beat-related features [58,59]. Often, the feature extraction process includes three 
stages: (i) Dividing each musical work into a set of short-term time windows, or frames, (ii) calculating 
a collection of features for each frame, and (iii) computing some sort of simple statistics (e.g., mean and 
variance) of each feature for all frames. The window size is chosen as a compromise between statistical 
significance and approximate stationarity of the data in each frame. Typical window sizes vary from 
20 to 100 ms [59,61,62]. 


3. Complexity Analysis and Visualization 


This section addresses the musical repertories of Frank Sinatra, Rolling Stones, Johnny Hallyday, 
and Julio Iglesias in the perspective of eight complexity indices and the MDS. The musicians were 
selected for their long and prolific careers, for representing different musical genres, and for singing in 
different languages. 

In a first phase, we apply Equations (1), (3), (4) and (6) to the time and frequency representations 
of the TS, ¥ and ||, respectively. Therefore, we characterize the musical works by means of the 
set of measures { Hr, PEr,Cr,CRr, Hs, PEs,Cs,CRs}, where the subscripts {T,S} denote the time 
and spectral complexity indices. For computing Hr and Hs the probabilities are obtained from the 
histograms of amplitudes of ¥ and ||, respectively, using 100 bins. For the PE, Cr, PEs, and Cs, 
we adopt the parameters d = 4 and t = 1, that were adjusted by means of numerical experiments. 
For computing CRr and CRs we adopt the Windows implementation of the gzip compressor, version 
1.3.12 (built upon the Lempel-Ziv coding algorithm L277). The variability of the individual quantities 
in the set {Hr, PEr, Cr, CRr, Hs, PEs, Cs, CRs} is analyzed and correlated with the artists’ musical 
careers. 

In a second phase, we consider that each individual index {Hr, PEr, Cr,CRr, Hs, PEs, Cs, CRs} 
captures distinct details of the musical works and that a more complete characterization is 
accomplished when using all indices simultaneously. However, since an 8-dimensional representation 
is not feasible, we adopt the MDS technique for dimensionality reduction and visualization. 
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3.1. Frank Sinatra 


Frank Sinatra (1915-1998) was one of the most popular singers of the 20th century. Sinatra’s 
musical style is close to ‘vocal jazz’, but there is still controversy and debate about this classification. 
In his artistic career of about 55 years, Sinatra recorded almost 60 studio albums and 300 singles, along 
with compilations and live albums. 

In this study we consider a total of 707 musical works included in 57 studio albums released in the 
period 1946-1993. The albums are ordered chronologically and referred to by the sequence i = 1,...,57. 
Therefore, we should note that the time lapse between two consecutive values of i is not precisely 
identical. Figure 2 depicts the evolution of the Hr and Hs (using the black marks + and 0, respectively) 
of the musical works versus the index of the album, i, where they are included. Given the dispersion 
of the Hy and Hg values, we group the musical works in windows of Ty = 5 albums centered at each i 
value (i.e., the window goes from i — 2 toi + 2), for improving the readability. Then, we calculate the 
25, 50, and 75 percentiles, and represent the results by means of three continuous lines. Numerical 
experiments showed that this width establishes a good compromise between limited volatility and 
accurate discrimination. Lower values of T,, increase the detail, but blur the charts, while higher 
values of T;, tend to filter too much the time details. We verify that there exist relationships between 
the evolution of Hr and Hs and the different periods of Sinatra’s artistic career, even knowing that 
these periods are neither rigidly defined nor absolutely consensual. For the other complexity indices, 
we reach similar results and, therefore, their representation is omitted here for the sake of parsimony. 


Album release year 
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Album index 


Figure 2. The complexity indices Hy and Hg versus i = 1,...,57 for Frank Sinatra’s music. 


The MDS is adopted for reducing dimensionality from an 8- to a 3-dimensional space, allowing 
a direct interpretation of the results. We start by constructing a 57 x 8 dimensional array, W = [wig], 
where wy, i = 1,...,57,k = 1,...,8, represents the median of the kth complexity index when grouping 
the musical works into windows of T, = 5 albums centered at each i value. Then, we calculate the 
dissimilarity matrices Ag = [6,4(uj,uj;)] and Ac = [dc(uj,u;)], {i,j} = 1,...,57, where d4 and 6c 
denote the arc cosine and Canberra distances between uj; = [wj : k = 1,...,8] and uj = [wy > k = 
1,...,8], respectively. The two distances are given by: 
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57 
yo UjKU jk 
6, (uj, uj) = arccos = , (10) 
57 5 57 5 
ye Mia) Wir 
k=1 k=1 


a Ujr| 
fo |Maxl + [eye 


6c (Ui, Uj) = (11) 

Other distances can be adopted, but several numerical experiments with distinct alternatives [63] 
confirmed that the arc cosine and the Canberra distances yield good results. Each of the matrices 
Aa and Ac is processed by means of the MDS for constructing the loci of objects that represent the 
evolution of complexity. 

Figure 3 depicts the MDS maps for Sinatra’s music, for g = 2 and q = 3, with Ay andi = 1,...,57. 
Figure 4a,b shows the corresponding MDS assessment charts. The Shepard diagram reveals a small 
scatter around the 45 degree line, demonstrating that there exists a good fit between the original 
and the reproduced distances. The stress plot shows that the maximum curvature of the line occurs 
close to gq = 2. Therefore, we conclude that q = 2 yields a good compromise between accuracy 
and readability of the locus of points, while q = 3 just leads to a marginal improvement, since the 
z-MDS coordinate carries reduced additional information. Alternatively, for taking advantage of 
present day computational visualization, we adopt a distinct 3-dimensional representation, with q = 2 
and the z coordinate of the map representing the albums’ sequence, i, interpolated with radial basis 
interpolation (RBI) [64] at each point with coordinates (x,y) produced by the MDS. The thin-plate 
spline @(€) = €? loge RBI function is considered, where the variable ¢ denotes the Euclidean distance 
between the points generated by the MDS for q = 2 and points in the xy MDS plane. Figure 5a,b 
depicts the results obtained for A, and Ac, respectively. The Shepard and stress diagrams are omitted 
here, since they are of the same type as the ones presented in Figure 4. 
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Figure 5. The MDS maps for Sinatra’s music. The (x,y) coordinates are generated by means of the 


MDS with q = 2 and the matrices: (a) Ay; (b) Ac. The z coordinate represents the albums’ sequence, 
i, interpolated with radial basis interpolation (RBI). 


We verify the emergence of five clusters, S;,r = 1,...,5. In the first, S; (ie., albums 1 <i < 12), 
the complexity varies strongly, meaning that the characteristics of the musical works evolved 
considerably. This cluster corresponds to albums released in the years 1946-1957. For S2 (13 < i < 25), 
the complexity has limited evolution and corresponds to albums recorded during the years 1957-1962. 
In the cluster S3 (26 < i < 34) the trajectory changes direction and has another large excursion, 
corresponding to musical albums recorded in the years 1962-1964. The cluster Ss (85 < i < 43) 
includes albums from 1964 up to 1967 and we verify that the complexity has a limited evolution. 
Finally, for the cluster Ss (44 < i < 57), another route occurs, smaller than the previous ones for S; and 
S3. Here, the complexity evolves slowly until the two last albums, consisting of duets (‘Duets I and 
II’), which explains the variation at the end of the career. It is also interesting to see that between two 
consecutive clusters S; and S;+4, (i = 2,3,4), we have always a trajectory tangle revealing the artist’s 
search for the new direction of work. 

We now analyze the musical repertory of Frank Sinatra by means of classical musical features, 
instead of general complexity indices. Therefore, each of Sinatra’s musical works (707 in total) is split 
into 50 ms non-overlapping frames, and a collection of 34 features is extracted for each time frame. 
After, for each feature, the average, the standard deviation, and the ratio between the average and the 
standard deviation are computed. Thus, each piece of music is characterized by a 34 x 3 dimensional 
vector, f;. Herein, we adopt the zero crossing rate, energy, energy entropy, spectral centroid, spectral 
spread, spectral entropy, spectral flux, spectral rolloff, mel frequency cepstral coefficients (13 values in 
total), chroma vector (12 values in total), and chroma deviation. For a detailed description about these 
features, interested readers can refer to [59]. It should be noted that a different set of features could 
have been used, since others are also available, and consequently many combinations are possible. 

We compute the 707 x 707 dimensional matrix A‘, = (54, (fi, f;)], where 6' (fj, fj) denotes the arc 
cosine distance between the feature vectors f; and fi, i,j =1,...,707. The matrix AY is used as the input 
to the MDS. Since the MDS technique outputs a large number of points, we post-process the results 
by (i) grouping the musical works into windows of T,, = 5 albums centered at each i value, and (ii) 
calculating the medians of the corresponding (x,y,z) MDS coordinates. Figure 6 depicts the resulting 
57-point 2- and 3-dimensional maps. Contrary to the previous experiments, in Figure 3, we do not 
see the emergence of any pattern. This means that in the perspective of this study, general complexity 
measures unravel characteristics somehow overlooked by specialized feature descriptors. While a 
systematic comparison of the two possible strategies, that is, the balancing between general indices 
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and specialized ones is of interest, hereafter we follow the first due to its superior performance in the 
present case. 
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(b) q =3. 
3.2. Rolling Stones 


The English rock band Rolling Stones was created in 1962. The original band included the vocalist 
Mick Jagger, the guitarists Keith Richards and Brian Jones, the bassist Bill Wyman, the drummer 
Charlie Watts, and the keyboardist Ian Stewart. Ian Stewart left the group in 1963 and Brian Jones in 
1969, being replaced by Mick Taylor, who remained until 1974. In 1975, the guitarist Ron Wood joined 
the band. The Rolling Stones quickly became the ‘bad-boys’ band, with an image of sex, drugs, and 
rebelliousness, in contrast to their contemporary band ‘The Beatles’. Their music was influenced by 
different styles from blues and jazz to dance and early rock-and-roll. The Rolling Stones are one of 
the most successful and acclaimed rock bands of all time. For more than 50 years, they released about 
30 studio albums along with several live albums and compilations. 

In the sequel we consider a total of 317 musical works included in 27 studio albums released in 
the period 1964-2005. 

Figure 7 depicts the 25, 50, and 75 percentiles of Hr and Hs versus i = 1,...,27, calculated as 
explained in the previous subsection. Again, we verify a relationship between the evolution of the 
indices and the different periods of the band’s artistic career. For the other complexity indices, we reach 
to similar results and, therefore, their representation is omitted here. 
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Figure 7. The complexity indices Hy and Hg versus i = 1,...,27 for the Rolling Stones’ music. 


Figure 8 depicts the 3-dimensional map of the Rolling Stones’ career. The (x,y) coordinates are 
obtained by the MDS with g = 2 and the dissimilarity matrices A, and Ac, while the z coordinate 
addresses the albums’ sequence, i, interpolated with RBI. The Shepard and stress diagrams are not 
represented, since they are of the same type as the ones presented in Figure 4. 

We verify that the complexity loci have two small tangles, S; (1 <i < 8) and S3 (18 <i < 22), 
corresponding to albums released between the years 1964 and 1966, and 1978 and 1986, respectively. 
These tangles intermediate two large excursions, S2 (9 < i < 17) and Sy (23 < i < 27), that include 
the albums released during the periods 1967-1976 and 1989-2005, respectively. In one hand, we can 
notice that the periods of complexity stagnation, S; and S3, comprise the early discography and the 
commercial success peak periods. On the other hand, the periods of strong complexity variation, 
S» and S4, include some troubled years and the entry of Ronnie Wood to the band, and the comeback 
and record-breaking tours that took place after the near break up. 
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sequence, i, interpolated with RBI. 


3.3. Johnny Hallyday 


Johnny Hallyday (1943-2017) was a French singer, songwriter, musician, and actor. He is 
considered the father of French rock and roll and sometimes he is referred to as the French Elvis 
Presley. Johnny’s artistic career lasted about 55 years and had plenty of musical success, especially 
in France and French-speaking countries. He recorded about 50 studio albums, as well as diverse 
compilations. He is well remembered for his spectacular live concerts with some shot of eccentricity. 

Figure 9 depicts the 25, 50, and 75 percentiles of Hy and Hs versus i = 1,...,34, calculated for a 
total of 325 musical works included in 34 studio albums released in the period 1961-2011. As mentioned 
for the previous artists, a relationship emerges between the complexity indices and the evolution of 


Johnny Hallyday’s career. 


Album release year 
-anmstmonraoanrmoernnonoowmrnanartnoerrtnwmaanwnonr DO 
© © GGOG&GooOoGBGRKRKEKRKRKRKRKR DD HHH HHA HGHHAH SSCS Se 
MPDRBDAHDADANHBA AAA HBAAAHAAAADAAAAHRA HAHAHAHA SSSSO 
a 3 3 

45 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
+ 
4- eye Sed 
toa , pesensaeeensseee ae ade: 
¢ re it 
4 : : é : goers—scege to 3 Be : # 
boone a a ae eared Seman St OE NS TE sa pS + 
5F FU a oe: t a * ee gat 5 “ ere 
eee eC ae ae eee eb toes yee 
+ = ate i ung eas aie 
a a ik 
+ 
3h : i e t —— Percentile 75 +} 
+ — Percentile 50 
~~ Percentile 25 
25 | | | | | | | | | | | | | | | | | | | | | | | | | | | | | I I I 
“- an metn OR DMOHT AMT MH OR DMRHOrTAMRYt HD OR DOROrA OM + 
ere ee Keer Ke KHON ANON N NO NNOUNOANANHOH OH HO OY 
-anmrtmonroanrmornnon oO owmraanartnoert nN waoanwnonrR DO 
© oOo GOoGGOHoHoEeRKRKEKRKRKRKR KR DD HHH DHA AHAB SESDCSSe 
DBD DADA HBA AHAAHAAAAAHAAHAAHAAAHAAAAAAAHAAHAHHAHSSSOSS 
a a a a a «CC 00 0 
44 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
8 ° 
° 8 9 
| 3 = 
4.2 6 @ 28 7 . ° 8 es 6 ee 
6__8 8 a 8 & ° 
pepe <SHse etl Mung eee ee ts Bae . 
4 3—g-4— + S SoreGiet ee escang Pd etka gerasecs Bio 
& 8 af ae 8 
pom | i eee eae er | 8 8B oS a $s 
30 8 o8 : eee oe cae 
3.84¢- ° 3} 3 2 8 tol o 
° ° 9 ° 
36 — $ «8 8 ~~ Percentile 75 |_| 
H —— Percentile 50 
~-— Percentile 25 
| | | | | | | | | | | | | | | | | | | | | | [ [ [ 
3.4 
rN MO tT OM OR DOHROrTr AMAT HM OR DWDAOrA OY - ao vt 
a | Oo 0 © 4 
Album index 


Figure 9. The complexity indices Hy and Hg versus i = 1,...,34 for Johnny Hallyday’s music. 
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Figure 10 represents the 3-dimensional map where the (x,y) coordinates are generated by the 
MDS with q = 2 and the dissimilarity matrices A, and Ac, and the z coordinate denotes the albums’ 
sequence, i, interpolated with RBI. 
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Figure 10. The MDS maps for Johnny Hallyday’s music. The (x, y) coordinates are generated by means 
of the MDS with q = 2 and the matrices: (a) Ay; (b) Ac. The z coordinate represents the albums’ 
sequence, i, interpolated with RBI. 


We observe the emergence of two main clusters. The first is a large tangle, S; (1 < i < 17), 
corresponding to albums released between the years 1961 and 1978. In this period, Hallyday recorded 
several French versions of American hits and French songs. The second, S2 (18 < i < 34), includes 
albums between the years 1981 and 2011. It begins in the early 1980s, when Johnny’s career seemed to 
be on the wane, and then evolves with a new breath triggered by the album “Rock’n’roll attitude”. 
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3.4. Julio Iglesias 


Julio Iglesias is a Spanish songwriter and singer. Iglesias’ career started in 1968 and has had plenty 
of commercial success and artistic recognition, with more than 300 million records sold, about 5000 
concerts for many millions of people, and dozens of awards worldwide. Iglesias is the most celebrated 
Latin music artist and one of the top 10 best-selling artists of all times. 

Herein, we consider a total of 629 musical works included in 59 albums released in the period 
1969-2007. 

Figure 11 depicts the 25, 50, and 75 percentiles of Hr and Hs versus i = 1,...,59, showing a 
relationship between the evolution of the indices and the different periods of the artist’s career. Such a 
relationship also emerges in other complexity indices, but their corresponding charts are omitted here 
for the sake of parsimony. 
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Figure 11. The complexity indices Hr and Hg versus i = 1,...,59 for Julio Iglesias’ music. 


Figure 12 depicts the 3-dimensional map generated with q = 2 and dissimilarity matrices Ay and 
Ac, with the z coordinate denoting the albums’ sequence, i, interpolated with RBI. 

We verify the emergence of four clusters. The first, S; (i.e., albums 1 <i < 22), coincides with 
the height of Iglesias’ success during the 1970s and 1980s of the twentieth century. In this period 
we observe that the complexity evolves as a tangle confined to a small region in the plane. For Sz 
(23 <i < 25), the complexity develops towards a new point, thanks to a few albums released in 1979. 
In the cluster S3 (26 < i < 34), the trajectory reaches another small tangle, corresponding to musical 
albums recorded between 1979 and 1982. This period precedes a fourth cluster, S4 (35 < i < 59), 
characterized by a large route and coinciding with albums released from 1983 up to 2007. Within 
this period, Iglesias started releasing many records tailored to suit American fans, including duets 
with some American stars. He then returned to the his Latin audience, including strengthening the 
relationship with his French followers, by releasing some French-language albums. In this period, 
Iglesias won the World Music Award and enjoyed major commercial success in Spain. 
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Figure 12. The MDS maps for Julio Iglesias’s music. The (x,y) coordinates are generated by means of 


the MDS with q = 2 and the matrices: (a) Ay; (b) Ac. The z coordinate represents the albums’ sequence, 
i, interpolated with RBI. 


4. Conclusions 


We adopted complexity, dimensionality-reduction, and visualization techniques for studying 
the music of several contemporary artists. The musical works were converted into digital format 
and represented in ‘mono’. The TS were assessed by means of eight distinct complexity indices. 
The 8-dimensional measurements were reduced to 2- and 3-dimensional by means of the MDS 
technique. The results revealed that the evolution of complexity is correlated with the artists’ musical 
careers. We conclude that the proposed indices represent reliable and assertive tools for assessing 
musical complexity. 
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